Local entanglement of multidimensional continuous- variable systems 
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We study the 'local entanglement' remaining after filtering operations corresponding to imperfect 
measurements performed by one or both parties, such that the parties can only determine whether 
or not the system is located in some region of space. The local entanglement in pure states of general 
bipartite multidimensional continuous-variable systems can be completely determined through sim- 
ple expressions. We apply our approach to semiclassical WKB systems, multi-dimensional harmonic 
oscillators, and a hydrogen atom as three examples. 
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I. INTRODUCTION 

It has been recognized quite recently that quantum 
entanglement is not just a profound feature of quantum 
mechanics but it is also a valuable physical resource, hke 
energy, with massive potential for technological applica- 
tions, such as quantum computation Jj, quantum cryp- 
tography [2] and quantum teleportation [3|, etc. How- 
ever, our understanding of entanglement is still far from 
complete despite current intense research activities. 

There are many reasons to focus on the entanglement 
of continuous- variable states [J, la, 0] > since the underly- 
ing degrees of freedom of physical systems carrying quan- 
tum information are frequently continuous, rather than 
discrete. Much of the effort has been concentrated on 
Gaussian states (i.e., states whose Wigner function is a 
Gaussian), since these are common (especially in quan- 
tum optics) as the ground or thermal states of optical 
modes. Within this framework, many interesting topics 
have been studied; for example, entanglement distillation 
for Gaussian states B, la, M, flOl|. multipartite entangled 
Gaussian states [ll|, E, [l3l ^^^ ent ang lement measures, 
such as entanglement of formation [IJ, ^M ^'^'^ logarith- 
mic negativity [16l [l7| . However one should remember 
that non-Gaussian states are also extremely important; 
this is especially so in condensed-phase systems, where 
harmonic behavior in any degree of freedom is likely to 
be only an approximation. Much less is known about the 
entanglement of these non-Gaussian states: while there 
is some progress in finding criteria for entanglement [l8| , 
there is little knowledge about how to quantify it. 

In two preceding papers [19l . |20|, we demonstrated 
how to use a specific type of projective filtering to char- 
acterize the distribution (particularly in configuration 
space) of entanglement in any smooth two-mode bipar- 
tite continuous-variable state. The approach is based 
on making an imperfect measurement of the 'position' 
of the system in configuration space, and then studying 



the entanglement remaining after the measurement. We 
showed how this approach could be used to map entangle- 
ment in different situations ,19.] , and that simple formulae 
exist for the entanglement in the limit where the region to 
which the system is confined after the measurement be- 
comes small (i.e., where the measurement becomes more 
and more accurate) :205]. 

In this paper we generalize these results to general (in- 
cluding multimode) smooth bipartite pure states. We 
first review the important results for two-mode states 
in 211 then generalize to multi-mode states in J jIIII Fi- 
nally in ijIVI we show examples of our approach applied 
to some systems in which analytical expressions for the 
energy eigenfunctions are easily obtained, before giving 
our conclusions in ^Vl 



II. TWO-MODE STATES 

We briefly recapitulate the definitions of essential 
terms and the known results for any smooth bipartite 
two-mode continuous-variable state. Let Alice and Bob 
share a state of two distinguishable one-dimensional par- 
ticles. Alice can measure only the position of her parti- 
cle (coordinate qA), Bob the position of his (coordinate 
qs)- They filter their state by determining whether or 
not the particles are found in particular regions of con- 
figuration space, and discard instances in which they are 
not. We refer to the resulting subensemble as the "dis- 
carding ensemble" . On the other hand if they choose not 
to discard the system when the particles are not in the 
desired regions, the resulting subensemble is called the 
"nondiscarding ensemble" . The entanglement £d in the 
discarding ensemble is related to the entanglement End 
in the nondiscarding ensemble by 
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where pat is the probability of finding Alice's particle 
within the region {qA '■ qA — o- l£ QA l£ QA + o-} and Bob's 
particle in the region {qs ■ qs — b < qs < qs + b}- We 
shall therefore focus on calculating £d, noting that £nd 
can be simply obtained from it; we show plots for both 
quantities for some of the systems discussed in mV\ 



A. Preliminary measurements on Alice's particle 
only 

If the initial state p is pure, so is po in the discarding 
ensemble. Suppose the initial filtering is performed only 
by Alice, by determining whether qA lies in the region 
QA — a < qA < 9yi + a, and all instances in which this 
is not the case are discarded. Now, since a is to be very 
small, Alice's original (before the measurement) reduced 
density matrix p^^{= Tr b[p]) in the neighborhood oi qA 
can be expanded (provided it is smooth in configuration 
space) as 

P^^Hqa; q'a) = P^^Kqa + x: qA + x') (2) 

= p(^) + pI^'x + p[fx' + p[fxx' + 0{x\x'\ 

where 



Within region —a < x < a, p\)' \s obtained by dividing 
equation ^ by the normalizing factor 2apQQ^ + 0(0"^). 

Now seek right eigenfunctions (/>„ of p^ within the 
allowed region: 



dx p'^^\x;x')(j)n{x) = A„0„(x). (4) 



Expanding 0„ as a power series 



(pnix) = an + bn{x) + —X + 0(a; ) 



(5) 
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(3) the eigcnfunction condition becomes (to order a^) 
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Expanding det(M — A/) to order a^ and equating to zero, 
we find two non-zero eigenvalues: 



From the eigcnfunction condition (j4|), wc find the third 
non-zero eigenvalue to be 



Ai = 



S^) 



3Poo 
A2 = 1-Ai. 



Mi' P'^ - P'^' P'^) 



(7) 



So to the lowest non-trivial order (a^), the eigenvalues, 

and hence the von Neumann entropy, of p]j are entirely 
determined by the quantity e = Ai. Specifically, the von 
Neumann entropy is 

5, - /!(£) ^ -[elog2(e) + (1 - e) log2(l - e)]. (8) 

To find the leading corrections to this result, we include 
all terms proportional to x^ or x'^ in the expansion ([2]) 
for p(^): 



SA) 



{x\x) 



(A) , (A) , (A) , 
Poo +/'10^ + /'01 ^ 



(9) 



^\{p'§x'^p[^x'^^2p^^^xx') 
+ \{p[fx^x' + pixx'^) 
^-^pifx^x'^+Oix^x'^). 



and then carry equation ([5]) to third order 

4>n{x) = a„ + bnX 



^CnX^ + ^dnX^ + 0{X^), (10) 



A., = 



{p'ofp^'p^^' 



y^Poo iPoi Pio "~ Pii Poo I 

+p'o^'pifp[^' + p'i'p'oi'pif 

-p[i'p'ofpii'~pii'p[i'pif)+Oia% (II) 



(A) (A) (A) 
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Therefore, the corrections due to higher eigenvalues, aris- 
ing from the higher-order terms in equation ([2]), affect e 
(and hence the entanglement) only to order a^. 



B. Preliminary measurements on both particles 

Now suppose both parties restrict their measurements: 
Alice's particle must lie in {qA '■ qA — a < Ia < qA + a|, 
and Bob's in {qB '■ qn — b < qs < qs + b}. In ^0| 
we attacked this problem by reducing it to an effective 
two-qubit one, for which exact results are available. How- 
ever this approach does not generalize so naturally to the 
multi-mode case, so we give here an alternative approach. 
From the argument above we know we can compute the 
entanglement from Alice's reduced density matrix p'^^' 
in the coordinate representation. Our first task, there- 
fore, is to evaluate this quantity once Bob has made the 
measurement of his particle. 



We do this by making a further Taylor expansion in- 
volving Bob's variables. We define 



Pnir 



QUI Qn2 Qns QUi 

p{qA,qB;q'A,(i'B) _ _ ■ 



(12) 



As we will see, to obtain the first nontrivial term in the 
solution we need all terms to first order in Alice's coor- 
dinates and to second order in Bob's: 

p{qA,qB;q'A,q'B) 

= p{qA + XA,qB + XB\qA + a:^, to + x'g) 

= Poooo + PimnXA + PowoXa + Pooioa^s + Poooia^s 

1 2 

+ 2(^00202^5^ +P00022^B ) 

+Piwoxax'a + Pwwxaxb + Pwqixax'b 
+Poiwx'aXb + Powix'ax'b + PoouXbx'b 
1 2 

+ -^{P1020XaXb^ + '^PlOllX AX bx'b + P1002Xax'b 

+P0120x'aXb^ + '2 PoillX A' xbx'b + Poio2x'ax'b ) 
+ OixA^x'/,XB^x'B''). (13) 

Alice's reduced density matrix is then found by writing 
rb 
'-b 



^{A) 



[XA 



1 /■'' 
y^A) = - dXBpixA,XB;x'^,XB) 

PJ-b 



26 

p 



[Poooo + xaPiooo + a^APoioo] 
+ ":r[Poo2o + 2pooii + P0002 

+ (P1020 + 2pioil -I- pia02)xA 
+ (P0120 + 2poill + P0W2)xa] 

+0{b'',XA^x'/). (14) 



where p is a normalization constant. By comparison with 
equation ^ and equating powers of Xi and j/i we can 
immediately identify the terms which appear in the ex- 
pression for e, and therefore determine the entanglement: 



(A) _ 26 
Poo — — [Poooo 



6 



(/O0020 + 2pooii + P0002)] + 0(6^ 



Pw^ = — [Piooo + -^(Pio2o + 2pioii + P1002)] + 0(6^) 
p 6 



(A) 

Poi 



.(A) _ 



26. 



[Poioo 



26. 



iPo 



120 



Pll = — [PllOO + -^(Pll20 

p 6 



2poiii+Poio2)] + 0(6^ 

2piiii+Piio2)] + 0(6^). 
(15) 



The leading (order 6^) terms in the numerator of the 
expression for e cancel — this is the reason why we need 



the density matrix to quadratic order in Bob's coordi- 
nates. The cancellation occurs because Alice and Bob 
(by hypothesis) share a pure state, and so 

p{qA,qB;q'AyqB) = iAqA.qBWiq'A^q'B) 



"3 



d" 



Qn 



dq'-^ 



?^"^ 



dqA"^^ dqB 



riq'A.qs) 



i^iqA,qB) 



(16) 



We can thus re-arrange the indices in a product of two 



terms as 



PabcdPef gh — PebgdPafch: 



SO in particular 



PiiooPoooo = PoiooPiooo- 



(17) 



(18) 



Hence the leading term in the numerator of e is of order 
6^, and the overall expression becomes 

„2u2 

[piioo(poo2o + 2pooii + P0002) 



loPoOOO 



+P0000(P1120 + 2piiii -t- P1102) 

— /Oiooo(poi2o + 2poiii + P0102) 

-P0100(pi020 + 2pioil + P1002)]. (19) 

Using equation ()17|) we can simplify this to obtain 
e 1 



a262 



^Poooo 



[pilOO/OOOll + POOOOPllll 



-pioooPoiii — PoiooPioiiJ 

1 



(20) 



l^Poooo 



[2piiooPooii + 2pooooPiiii — PioooPo 



111 



—poiooPioii ~ PooioPiioi — PoooiPiiio]- (21) 

The first form (|20p is slightly more compact, while the 
second form (pij) makes it clear that the coordinates of 
Alice's and Bob's subsystems are treated equivalently, 
as required. The von Neumann entropy, and hence the 
entanglement (since this is still a pure state), is then 
Sv = h{e) as before. 

We know, from the arguments leading to equation (jlip . 
that the leading correction to this result is 0(a*), and we 
should expect from the symmetry between Alice's and 
Bob's systems that it is also 0(6"*). We have explicitly 
computed the correction and this is indeed the case: the 
result is given in Appendix |^ The third eigenvalue A3 
measures the extent of the breakdown of our approach. 
We note that it is of order a'* 6^, and therefore does not 
affect the expression of e, which is of order a^b"^ . 



III. MULTI-DIMENSIONAL SYSTEMS 

A. General approach 

Consider first the case in which only Alice makes 
preliminary measurements. If Alice's system is two- 
dimensional and she localizes the particle so —ai < Xi < 



+ai, i £ {1,2}, one can find the eigenvalues of p*^"^^ by a 
straightforward generahzation of the methods in JjII Al 
Once again we find that there are only two non-zero 
eigenvalues to order af: 



Ai = 



A, 




.(A) d^p^^^ dp^^^ dp' 



{A) 



dqAAdq'j^i dqA,i dq'j^^ 



(22) 



where i goes over the two spatial dimensions of Al- 
ice's subsystems, H.T. stands for higher-order terms and 

We now argue that this property holds irrespective of 
the dimensionality of Alice's system, as follows. The 
entanglement must be invariant under exchange of the 
axis labels, and under all transformations of the form 
ai — > —ai. The only possibilities consistent with these 
requirements are 

Ai = l-^t.a2; \2 = Y,Uah A3,A4... = 0, (23) 



only the first form (|23p is consistent with equation (|22p . 
Therefore, the form of the non-zero eigenvalues must be 



^1 -i:^(^2[p 

+H.T. 
A2 = 1-Ai 
A3 = + H.T. 



dqAadq'j^.^ dqAa Na^i ) 

(25) 



where i now goes over all the dimensions of Alice's sub- 
systems. 
Define 



P(i,j;nin2n3ni) 



QUI Qn-2 



Qns 



Qrii 



dqA.'^' dq'A^ ^tej"' ^9s,, 

P{<^AMA.^B,q^B)\q^^^q^y (26) 



where i (j) represents one of available dimensions of Al- 
ice's (Bob's) subsystem. If the state /3(qA,q^,qs,qs) is 
pure, we have the following relation 



Ai = l-^i,a^; A2 = iia?; A3 = ^303, . . . , (24) 

i 

where the ti are arbitrary constants. Furthermore the 
eigenvalues must reduce to the known forms for one- and 
two-dimensional systems if all other a^ are set to zero. If 
we keep ai and 03 non-zero, sending all others to zero. 



P(i,j;nin2n3n4)P(ij";n5n6"7"s) ~ ^(i J';"5"2n7n4)^(« J':nin6n3"s) ' 

(27) 
From the previous analysis that led to equation pS]) 
for a pure two-mode state, we know we can extend equa- 
tion (I25p to a pure multi-dimensional bipartite state 
p(q^,q^,qs,q'^) for the case where both parties make 
preliminary measurements on their particles by making 
the following substitutions: 



J 



(-4) 

Poo = 



(-4) 
PlO 



E 
E 



p'.f - E 



p[f = E 



l^n, 


, 2b,, 


v 


P 


fn, 


, 2b,, 


v 


P 


m, 


, 2b,, 


v 


P 


(n, 


,2by 



Piij;0000) + -^{P{ij -,0020) + 2/?(ij;0011) + P(ii;0002)) 

P(ii;1000) + "^(P(ii;1020) + 2/?(ij;1011) + P(ij;1002)) 

&2 

P{ij;OWO) + "^(P(ii;0120) + 2p(ij;0111) + P(ij;0102)) 

P(ii;1100) + "^(P(ii;1120) + '^P{ij;llll) + P(ii;1102)) 



H.T. 



H.T. 



H.T. 



H.T. 



(28) 



where j and j' go over all the dimensions of Bob's sub- 
system and p is an appropriate normalization constant. 

Therefore, to the lowest order in a and b, Ai in equation 
(l25t becomes 



Ai = E 



«?62 



»,J 



18p(ij;oooo)' 



"|P(ij;1100)[P(ij;0020) 



(29) 



+ 2p(ij;0011) + P(ij;0002)] 

+P(ij;0000)[P(ij;1120) + "^PdJ;!!!!) + P(i,i:1102)] 
— P(ij;1000)[P(ij;0120) + '^P(i,j;0111) + P(i,i:0102)] 
~P(ij;0100)[P(ij;1020) + 2p(,;j;ioil) + P(i ,i;1002)] | ■ 

This can be further simplified by using equation (|27p to 



obtain 
Ai = 



E 



^,3 



9p(i,j;0000) 



r[P(i,i;1100)P(i,i;0011) 



+P(i,i;0000)P(i,i;llll) 
~P(i,j";0100)P(i,i;1011) 



P(i,i;1000)P(i,i;0111) 



(30) 



Again the entanglement is completely determined by 
Sy = h{e), where e = Ai as before. 



B. Concurrence and negativity for general 
bipartite multi-mode pure states 

In a similar way, we can generalize our previous expres- 
sions [2(| for the concurrence 21 1 and negativity [22, _23J 
of the system after the preliminary measurement has 
been made. 

For an ni ® n2 (jii < ^2) bipartite system, where ni 
and n2 are Hilbert space dimension for two subsystems 
respectively, the generalized concurrence of a pure quan- 
tum state V is defined by [2J| 



C'(|^))=4^ A,„A„, 



(31) 



where \/Xm {m — 1, . . . ,ni) are the eigenvalues of the 
reduced density matrices p^"^^ and p^^\ Additionally, 
the trace norm of the partial transposed density matrix 
with respect to Alice's subsystem turns out to be 



nTA\ 



(Ev^) 



(32) 



From this we can determine the negativity, which is de- 
fined as 



mp)^ 



o^^W-l 



(33) 



As we argued earlier, the reduced density matrix in the 
discarding ensemble has only two non-zero eigenvalues 
(Ai and (1 — Ai)) to the lowest order so we then have 
from equation (PT|) : 



4 2_^ \n\ 



4Ai +H.T. 

((E v^)' - 1) 



(34) 



where we have used X^m Am = 1- Therefore, we have 
proved that in the limit of small ai and hj , for any multi- 
mode bipartite pure state t/i. 



C{iP) = 2M{ip) = 2\/i. 
Specifically, the squared concurrence is 



c'-E 



2aihj 



"J 

/ .'^iy 



V- 



av 



9V d^ 



dqA,idqBj dqA,i dqs, 



(35) 



(36) 



where i goes over all dimensions of Alice's subsystem and 
j of Bob's subsystem. C? is the squared concurrence 
associated with the degrees of freedom i and j. Note 
that Cij oc ttibj, consistent with the existence of a well- 
defined local concurrence density for two-mode systems 

Note also that the concurrence is made particularly 
simple by writing 



in which case 



C2 



E 



V' 



4a^ 
9 



-,-s 



d^S 



dqA,idqB,j 



(37) 



(38) 



From this, we see that if S is quadratic in the coordinates 
(i.e., the state is a Gaussian) the local entanglement is 
constant; on the other hand whenever S' is a linear func- 
tion of the coordinates, the local entanglement is zero. 



Nodes in the wavefunction 



Evidently S in equation (|37|) diverges near nodes of 
the wavefunction, so that for a fixed ai and bj the con- 
currence given by equation ([55]) also diverges (like 1/|V'P 
as IV'I -^ 0). It is important to realize that this diverg- 
ing quantity refers to the entanglement in the discarding 
ensemble (i.e., in the sub-ensemble conditional on find- 
ing the particles in the chosen measurement region — see 
equation [T]) , and that even in this ensemble our expres- 
sion applies only in the limit of very small measurement 
regions. We now show that the discarding entanglement 
always remains finite provided we keep within the domain 
of validity of our approach. 

The extent of the domain of validity follows inevitably 
from our Taylor-series approximations for the wavefunc- 
tions (or density operators — see equation (|16[) ). which are 
valid only close to the chosen reference point (qyi,qB)- 
The requirement that the second term in this expansion 
be small compared with the first is 



dip 

^qA.i 



tti < V(qA,qB) 



a, < 



•0(qA,qB) 

dtp/dq 



AA 



(39) 



and similarly for bj] therefore, the domain of validity 
shrinks to zero near a node in ip. Equivalently, if this 
condition is not satisfied it leads to the breakdown of the 
isomorphism of each mode to one qubit described in ,20| . 
One way to understand the behavior of the entangle- 
ment near points where the wavefunction vanishes is to 
satisfy equation p9p by writing the maximum valid re- 
gion size as 



a,- = a- 



dtp/dq/ 



(40) 



where a <C 1 is a small parameter, and similarly for 
(We assume here that the derivatives are not also 



lMAX 

j 



zero near the nodes.) We further define three quantities 
ki, kj, and kij by 



and similarly 



dqA,idqB,j 
dtp 

dqA,i 

dtp 



= kitp; 
== kjtp, 



(41) 



so af^^^ki = hf^^kj = a. From equation dSS]), if we 
choose ai = a^^-^, hj = h^^^ near a node where kikj ^ 
fcy, the expression for e reduces to 

4 

EMAX^^y. (42) 

Therefore e (and hence also the localized concurrence and 
entanglement) is cut off near the node at a finite value 
that depends on the choice of a. 






dp dp 
dQ, dQ', 



(46) 



Therefore, equation ()25p is invariant under the gener- 
alized orthogonal transformation (H^)) . It follows that 
equation ((30)l . and hence the local entanglement, are also 
invariant under these local transformations. 



2. Non-local transformations 

We now consider some transformations which mix Al- 
ice's and Bob's coordinates — specifically, those that make 
the system separable. That is to say we look for a new 
set of coordinates 



Xk — 2_^ TikXi 



(47) 



D. Transformation of coordinates 

We now discuss the behavior of our expressions for the 
local entanglement under various coordinate transforma- 
tions. 



1. Invariance under local transformations 

We would expect that the definitions of our local en- 
tanglement would remain unchanged if wc made a local 
redefinition of our coordinate axes (possibly accompa- 
nied by changes in the measurement region) . To see that 
this is the case, consider the following transformation of 
Alice's coordinates: 






(43) 



where O is an orthogonal matrix {OO^ — 1) and the 
sum goes only over the other coordinates of Alice's par- 
ticle. Ai is to determine the length of the measurement 
region for new variable Qi. Note that if aj = Ai = a^i, j 
(i.e. both measurement volumes are hypercubes with the 
same dimensions) then (|43l) reduces to a simple orthogo- 
nal transformation of Alice's coordinates. 
Now 

d ^.dQ. d ^^ . d ^^^^ 



dqr 
We then have 



y^A-^— 



dqi dQj 



3 ■' 



E 



2 d^P 



T = ^Oifi.kAjAk 



d'^p 



dQ.dQ' 



E^: 



2 d^P 



(45) 



such that the wavefunction factorizes as 

tP = Y[MXk). 



(48) 



Note that the sum over i in ([TT]) runs over all coordinates 
of the system (both Alice's and Bob's). In this situation 
it does not make sense to consider any accompanying 
change in the shape or size of the measurement region, 
which we continue to define in terms of the original co- 
ordinates and to describe by {at} and {bj}. 
Therefore, 



d'^tP 
dxidxj 



kk' 



9V 



dXkdXk' 

ip d^tpk 
^k dXl 



(49) 



E T^kT, 



3k' 



k^k' 



tp dtpk dtpk' 
tpktp'^. dXk dXk' 



and similarly 



dtp_dtP_ _ v^„ ^ 



kk' 



It follows from equation 



E 



{a.bjY 



k ^^ 



tp'^ dtpk dtpk' 
tpktp'^ dXk dXk' ' 

that 

d^i^k l_ f dtPk 
i^k \dXk 



(50) 



dXi 



(51) 
where the second term inside the modulus signs comes 
from the part of ([50]) having k = k' . In terms of 
the logarithms of the separable wavefunctions 5*^ = 
-log[V'fe(^fe)]), we have 



= E 



(a^bj)^ 



J-lk-Ljk f^Y2 



dxl 



(52) 



One important special case of this result is the trans- 
formation to normal coordinates in a harmonic system: 
if the potential can be quadratically expanded about an 
energy minimum, the transformation to normal coordi- 
nates takes the form of equation (jlT]) with 



(53) 



where O is an orthogonal matrix. 

3. Relative coordinates 



A closely related example is the transformation to 
center-of-mass and relative coordinates. (Here we as- 
sume that the particles live in the same physical space, 
and hence that the dimensions Na and Nb are equal.) 
If Alice's particle and Bob's particle have masses uia 
and niB respectively, we define r^ = qf — qf and Ri = 
{fi/mB)qf + {n/mA)qf where fj. = mAmB/{mA + tub) 
is the reduced mass and i goes over all dimensions of the 
system ({x,y,2:} in three-dimensional system, for exam- 
pie). 



E 



O-ibj 



M 



(9V 



d 



niB dRi 
d 



-)( 



A* 



d 



dip dtp 



dr. 



niA dRj 



d 



TUB dRi dri niA dRj drj 



)V; 



(54) 



where i and j run over all the dimensions of the system. 
In many cases, including most importantly the case 
where there is no external potential, the wave function 
'(/'(R, r) can be decoupled into a center-of-mass part xC^) 
and a relative- motion part f{r): 



7/;(R,r)=x(R)^(r). 



If we write 



(p(r) 



- p-s^W 



X 



(R)=e-^x(^) 



(55) 



(56) 



then the entanglement takes the particularly simple form 



E 



alb] 



d^S^jr) , fi^ 



d^S. 



dvidrj mAmB dRidRj 



(57) 



For example, if x(R) is a free-particle plane wave x(R) = 
gikoR^ its contribution to the entanglement £d is zero; if 
x(R) is a Gaussian wave packet with wave number ko 
and real-space width Rq: 



V^(R,r) = (4j)V4e-RV«;^e<koR^(r), 



^nRl 



the expression for e becomes 



2l2 



afb 



9 



d^S^ir) 



2fi^ 



drjdri 



mAmB-n^ 



Rr' 



(58) 



(59) 



V{r) 




FIG. 1: Diagram of a potential well illustrating the different 
regions discussed in the text. 



IV. EXAMPLES 

In this section we apply our method to some easily 
soluble examples: first to wavefunctions that (while re- 
maining pure states) are semiclassical in the sense that 
the potential varies slowly on the scale of the de Broglie 
wavelength, so WKB methods are applicable, then to en- 
ergy eigenstates of harmonically-interacting particles in 
arbitrary dimensionality, and finally to bound states of 
an electron and proton (i.e., to the hydrogen atom). 



A. The semiclassical case: one-dimensional WKB 
\vavefunctions 

Consider two particles moving in one dimension with 
an interaction potential V{r) that depends only on the 
relative coordinate. Neglecting center-of-mass contribu- 
tions, the entanglement can then be calculated from the 
relative wavefunction ^p{r) using equation (|57p . If V{r) 
is a slowing varying function of r, we can use the WKB 
method to find ip{r). 

We consider an interaction with a single potential well 
(shown schematically in Figure [1]), so the system moving 
in a bound state with energy E has just two classical 
turning points. For the classically allowed region with 
E > V (region 2 of Figure [1]) , the classical momentum 
at r is p{r) = ^j2m{E — V{r)) and the corresponding 
wavefunction can be expressed as 



^,— (r) 



2A 



VpF) 

so that the local concurrence is 
ab 



p{r') dr' + 



Co = 



ih'p{rf 



2csc^ 



1 r^^ 

h 



p{r')dr' + - 



^1 I ^ d^PJr) 



f dp{r 



\ dr 



ri < r < r2 
(60) 

p{r)' 
(61) 



+2acot 



-/ p(/)d/ + - 



V{r) 



2 dv(r) 



dr 



The oscillatory structure of the wavefunction, arising 
from the interference between right- and left-moving 
travelling waves, produces nodes at which the entangle- 
ment in the discarding ensemble for fixed a and b diverges 
(but remains finite provided we remain within the domain 
of validity of (l61])--see ^Q. 

Note also that the entanglement contribution from the 
first term in (|6ip is non-zero even where V{r) (and hence 
p{r)) is constant. 

For E < V (region 1 and region 3 of Figure [1]), 
we express the wavefunction in terms of the local 
momentum on the inverted potential surface p{r) — 
yj2m{V{r) — E). The wavefunctions are respectively 



^r^^(r) 



^f^^(r) 



(-1)"A 



exp 



A 



VW)\ 



exp 



\p{r')\dr' 



|p(r')|dr' 



r < ri; 
(62) 
r > r2, 
(63) 



where n is the number of nodes in Region 2. Correspond- 
ingly, the concurrences are 



Ci 



~ab 



m\p{r)\' 



2\p{r)\ 



d\p{r)\ 



C. = 



ah 



3n\p{r)\' 



d\pir)\ 
dr 

2\p{r)\ 



d\p{r)\ 
dr 



dr 




h\pir)\ 


d' \p{r)\ 


dr^ 


d\p{r)\ 




dr 




n\p{r)\ 


d' \p(r)\ 



(64) 



(65) 



Note that in this case (by contrast to the behavior in re- 
gion 2) if there is no force, p{r) is constant, and hence 
there is no entanglement. It is interesting that the bound- 
aries between these different behaviors of the entangle- 
ment correspond to the classical turning points. 



B. Multi-dimensional harmonic oscillators 

Consider first a system of two one-dimensional har- 
monic oscillators of masses m^ and ms , having identical 
frequencies u>, and coupled by a spring constant K; the 
Hamiltonian is 

H^HA + HB + \K{XA + XBf. (66) 

Transforming to center-of-mass and relative coordinates, 
the eigenstates are simply 



v/^/^2"«2"'-r^fl,!r^r!i?oro 
R 



,-RV2RI 



e-^ '^^^Hr.A-fr)HnA-), (67) 
-Ko rQ 

where nn and n^ label the excitations of each coordinate, 

i?o = ^n/{Mu), ro = ^[h/iil^/^TWf^), and H,,{x) 
is the Hermite polynomial. 

If Alice and Bob each possess an oscillator, the en- 
tanglement between their subsystems given by h{t) can 
be determined from equation (|54p : for example, for the 
ground state: 

a%'^{mAmBrl - M'^Rlf 



^M^r^Rl 



IK 



QM^h^ 



{mArriBUJ — Mu\ / \- uj 

V f^ 



2^2 



(68) 



where M — niA + tub- Note that the ground state is 
Gaussian, so e is constant, as expected. 

In Fig. [21 we plot the probability distributions and 
entanglement £ (in the discarding ensemble — center col- 
umn, and nondiscarding ensemble — right column) for the 
ground state and some excited states. Note that the 
ground state (a) is a Gaussian state so the discarding 
entanglement is constant and the left and right plots are 
proportional to one another; this is no longer true for 
the other (non-Gaussian) states, for which there are also 
nodes in the wavefunctions. We therefore show the en- 
tanglement in both ensembles cut off at the maximum 
value determined by equation (|42l) . 

For general multi-dimensional oscillators, the wave- 
function becomes a product over the normal modes Xk of 
one-dimensional harmonic oscillator wavefunctions. The 
entanglement is determined by these normal-mode wave- 
functions through equation ([52]) . (Note that in the one- 
dimensional example considered above, the normal coor- 
dinates are the same as the relative and center-of-mass 
coordinates.) 



C. The hydrogen atom 

We next consider the entanglement between the elec- 
tron ('Alice's particle') and the proton ('Bob's particle') 
in a hydrogen atom. For simplicity, the sizes of the mea- 
sured regions are assumed to be the same for all dimen- 
sions {x, y, 2;}, i.e. ai = a and bi = b. First, consider the 
case where there is no center-of-mass motion. Instead of 
directly applying equation (|57p . we transform the coor- 
dinates and the equation into to spherical coordinates: 



d 

drx 


. ^ j^d ^ cos0cos</) d CSC 6* sin d 
dr r do r dcp 


d 
dry 


d cos 6* sin d csc6'cos(/> d 

= smOsmcp—- -i -— H — 7 

dr r do r dcj) 


d 
drz 


-<~w (-' 





(A) Vo,o(-R,r) 




(C)^i,3(i?,r) 

FIG. 2: Probability density (left plot), local entanglement £d 
in the discarding ensemble (center plot) and local entangle- 
ment End in the nondiscarding ensemble (right plot) for three 
pure states of the two-oscillator system: (A) ur = 0, Ur — 0; 
(B) Ur = 1, Ur = 1; (C) Ur = 1, Ur = 3. The characteristic 
lengths of the problem are ro = 2 and _Ro = 4 in all plots, 
and all plots are for a = 6 = 0.1. The cut-off points for plots 
of £d and End are determined from emax in equation (|42p 
with a 
defined in equation 



0.1; specifically, £^^^ 



P^b^^H^MAx), where h is 



The ground state is 

^ioo(r,0,0) = (^)i/2e-'^/'^°, 
where ao is the Bohr radius. In this case, 



2( 



3anr 



(70) 



(71) 



Interestingly, this expression indicates that the entangle- 
ment En for the ground state of a hydrogen atom falls off 
with distance in exactly the same way as the electrostatic 
force between the electron and the nucleus. 

If we include a center-of-mass part to the wave function 
with a Gaussian form as in equation ([55]) . we obtain 



2a^^ 



Rt{mj 



niB) 



9i?;^a2(TO^ + mB)4r2 

— 4i?Qaom.A'7i_B(m^ -I- mB)^r + Qa^rn^m^r^ 



(72) 




FIG. 3: Probability density (left plot), local entanglement Ed 
in the discarding ensemble (center plot) and local entangle- 
ment End in the nondiscarding ensemble (right plot) for the 
relative wavefunction ip2\o(r,6,(j)) of a hydrogen atom. All 
plots are for a = b = 0.1. The cut-off points for plots of Ed 
and End are determined from emax in equation (|42p with 
a = 0.1; specifically, E^j^'^ = p^ '^h{eMAx), where h is de- 
fined in equation ((SJ. 



The first term is the component noted previously, decay- 
ing in the same way as the atom's internal electrostatic 
force; in addition there are two new contributions from 
the localization of the free-particle wave function. Of 
these the third term corresponds to the spatially con- 
stant entanglement of the gaussian center-of-mass state. 
Excited states of the atom can also be analyzed, by 
substituting the most general form of the relative wave 
function (fnim ir^ ^j 0) of a hydrogen atom into equation 
(|57l) after it has been transformed to spherical coordi- 
nates. The excited states have nodes in the wavefunction, 
which have to be treated as discussed earlier. We show 
the corresponding probability distribution, and entangle- 
ment £ (in the discarding and nondiscarding ensembles) 
in Fig. El 



V. DISCUSSION AND CONCLUSIONS 

Our approach allows us to analyze the distribution 
of entanglement after imperfect local position measure- 
ments in any smooth bipartite pure state. Equations 
(PS)) and (|38p are our main results, allowing us to cal- 
culate the concurrence in terms of simple derivatives of 
the wavefunction. Equation (|5T|l allows us to express the 
entanglement in the same local region in terms of an ar- 
bitrary linear transformation of the coordinates, and (|57p 
treats the important case where the motion separates into 
center-of-mass and relative coordinates. 

The three examples of exactly integrable systems that 
we have discussed show a number of common features. 
First, there is generic behavior near nodes in the wave- 
function. There is an apparent divergence in the entan- 
glement in the discarding ensemble for a fixed region size, 
but this does not mean that large amounts of entan- 
glement can be extracted from the continuous-variable 
wavefunction once the system has been localized in this 
region. Our expressions for entanglement are always true 
only in the limit of small region sizes, and their domain 
of validity shrinks as we approach a node; the discarding 
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entanglement remains finite so long as we take care al- 
ways to remain within this domain. Furthermore, when 
we measure the locations of the particles we are unlikely 
to find them near a node in the wavefunction, so the 
probability factor in equation ([T]) further suppresses the 
non-discarding entanglement relative to the discarding 
entanglement. 

As the size of the measurement regions increases, our 
approach starts to break down because more than two 
eigenvalues of the reduced density matrix become impor- 
tant. We have explicitly computed the extent of this 
breakdown, giving the lowest-order corrections to our 
main results in Appendix [X] 

As pointed out in §111 D 3[ free-particle wavefunctions 
do not give rise to any local entanglement. We have 
shown how our entanglement expressions are transformed 
when moving to other coordinates (e.g. center-of-mass 
and relative coordinates); however, it is important to re- 
alize that the entanglement we quantify is still between 
the original subsystems. The transformation is only done 
for the convenience of the calculations. 

Our results for the WKB wavefunctions and for the 
hydrogen atom suggest an intriguing link between the 
interaction force and the local entanglement, but the ex- 
act details of the relationship and its generality need to 
be further explored. We also note that in this paper we 
have only considered pure states; the application of our 
approach to the mixed states will be discussed in another 
paper. 



APPENDIX A: CORRECTIONS TO THE LOCAL 

ENTANGLEMENT AFTER TWO-PARTY 

PRELIMINARY MEASUREMENTS 

The third eigenvalue of Alice's reduced density ma- 
trix in the discarding ensemble when both parties make 
preliminary measurements can be found by making the 
following additional substitutions in equation (|11[) : 



(A) 
P20 



P 



[P2000 + 7r(P2020 + 2p20ii + P2002)] + 0(6^); 



(A) _ 
P02 — -[P02OO 



6 



7r(po22o + 2po2ii + P0202)] + 0(6^); 




P21 = -[P2ioo + 7r(/02i20 + 2p2iii +P2102)] +0(5^); 
p 



6 



7 T 2 

P22 = ^[P2200 + -p(P2220 + 2p2211 + P2202)] + 0{b^). 

(Al) 



This gives 



where the denominator is 



\de ' 
^3 



(A2) 



\ de 



120(pooo2PoiooPiooo + 2pooiiPoiooPiooo 
+P0020P0100P1000 + P0000P0102P1000 
+2pooooPoiiiPiooo + P0000P0120P1000 
+P0000P0100P1002 + 2/3ooooPoiooPioii 

+P0OOOPOIOOP102O ^ 2/3oOOOP0002PllOO 
— 4pooooPooiiPiioo — 2poOOOP0020PllOO 
— P0000P0000P1102 — 2pooooPooooPiiii 
— P0000P0000P1120), (A3) 



and the numerator is 



\ nu 
-^3 



77(P021lPll20P2002 - 
54 

— POIIIP1220P2OO2 
+PO2O2PII20P2OII 

— POIO2P1220P2OII 

— PO220PIOIIP2IO2 
+POO20P12IIP2IO2 
~P0220Pl002P2111 
+POO20P12O2P2III 

— PO2IIPIOO2P2I2O 
+PO2IIPIOOOP2I22 
+POI20PIOIIP22O2 
~P0020PllllP2202 
+POI20P1OO2P22II 

— POO20PIIO2P22II 
+POIIIPIOO2P222O 
~P001lPll02P2220 
~P011lPl000P2222 
+POOIIPIIOOP2222 



POI20P12IIP2OO2 
" P0220PIIIIP2OO2 

- POI20P12O2P2OII 
' PO220PIIO2P2OII 

- P02IIPIO20P2IO2 
" POOIIP1220P2IO2 
' P02O2PIO20P2III 
" POOO2P1220P2III 

- POOO2P12IIP2I2O 
' POOOOP12IIP2I22 

- POIIIPIO20P22O2 

- POOIIPII20P22O2 

- POIO2PIO20P22II 
' POOO2PII20P22II 

- POIO2PIOIIP222O 
' POOO2PIIIIP222O 
' POIOOPIOIIP2222 

" P0000PllllP2222)- (A4) 



P12 — -[P12OO 



7r(pi220 + 2pi2ii + P1202)] + 0(6^); 
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